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O ' I. INTRODUCTION 

<N : 

Since the pioneering work on coherent averaging effects by Haeberlen and Waugh Q, the use of tailored pulse 
sequences for manipulating the effective Hamiltonian experienced by a target quantum system has developed a solid 
tradition in nuclear magnetic resonance (NMR) In particular, within the context of NMR quantum information 

processing (QIP), decoupling and refocusing techniques provide the basic tools for enforcing Hamiltonian evolutions 
| that correspond to quantum logic gates between selected spins The principles underlying these techniques, along 
with the powerful formalism offered by average Hamiltonian theory (AHT) jl], have been recently extended beyond 
t-H ' the NMR domain, and suggestive applications have resulted in various directions within QIP. On one hand, ideas 
from NMR decoupling motivated a "bang-bang" S control-theoretic framework for generic open quantum systems 
H , which paved the way for the development of quantum symmetrization procedures and quantum error suppression 
strategics for QIP |^|-|l0|. On the other hand, the application of active dynamical control in the bang-bang limit 
proved a valuable tool for engineering the evolution of coupled quantum subsystems leading to various schemes 
_^ | for universal simulations of both closed-system p"T[ and open-system dynamics p"2| ] . 

-H ■ So far, in spite of the pervasive role played by quantum coding in QIP, the application of active refocusing and 
decoupling methods has been primarily thought of in terms of the basic physical degrees of freedom. Two exceptions 
are a proposal by Wu and Lidar for applying recoupling schemes on encoded qubits governed by exchange-type 
Hamiltonians pUf , and an implementation by Fortunato et al. of encoded refocusing to gain universal control on 
a decoherence-free qubit JTij . It is the purpose of this note to further comment on the significance of dynamical 
control methods as directly represented in terms of encoded degrees of freedom, by continuing the investigation of the 
interplay between quantum coding and decoupling techniques undertaken in || , and expanding the basic arguments 
sketched in fl4f| . While decoupling methods have been already shown to enable, in principle, to synthesize effective 
evolutions supporting noise-protected, encoded degrees of freedom ||,[l5|], a different perspective is taken here - by 
imagining a pre-selected encoded structure, and by looking at the evolutions that can be enforced through encoded 
decoupling sequences. 

I envision two prospective types of applications of encoded dynamical decoupling methods within QIP. Similar to 
their un-encoded counterparts, these include control of both closed-system Hamiltonian evolutions and open-system 
noisy evolutions - all evolutions being however restricted to an underlying coding space that is preserved throughout. 
After summarizing the basic facts about decoupling in Sect. II, encoded dynamical decoupling is introduced in Sect. 
III. I then discuss the two relevant areas of application based on simple representative examples in Sect. IV and V, 
respectively. A brief summary concludes in Sect. VI. 



II. BASICS OF DYNAMICAL DECOUPLING 



A bang-bang (b.b.) quantum control problem is concerned with characterizing the effective evolutions that can be 
engineered by repeatedly interspersing, according to various possible schemes, the natural dynamics of a quantum 
system with full-power, instantaneous control operations (b.b. controls) Let S denote the target control 

system, defined on a (finite-dimensional) state space Tts, dim(Ti.s) = N, N — 2™ for qubit systems. In a general 
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open-system setting, S is coupled to an environment E via an interaction Hamiltonian Hse- The control problem 
can then be formulated in terms of the following data: 

• H, the natural Hamiltonian of the joint system, H = Hs <E> 1b + I5 <£> He + Hse, determining the free unitary 
evolution Uo(t) — exp(-iHt) on the joint state space Hs 8> He- One can write Hse — Yl a ^ a ® ^a, ^ or appropriate 
system (S a ) and environment (E a ) operators. The linear subspace Af = span{5 Q } of noise-inducing couplings is 
referred to as error space. Without loss of generality, both Hs and the error generators in AT can be assumed to be 
traceless. Detailed knowledge of Hs or Af may or may be not explicitly available from the start. 

• Gb.b., the set of all b.b. operations that can be effected on S. As it is conceivable that P^ 1 = € Gb.b. if P € Gb.b., 
Q can be associated with a subgroup of the group U (H) of unitary transformations on the system alone; 

• G, the discrete subset of b.b. decoupling operations, G — {Uk} Q Gb.b., with k <E K for some finite set of indexes 
with order \JC\ = \G\; 

• T c , the relevant control time scale (cycle time), associated with the duration of a single control cycle; 

• {rfe}, the set of relative temporal separations between consecutive b.b. operations, = Atk/T c > 0, in terms of 
the free evolution intervals Atfc, and X)fce/c Tk = T 

AHT provides a general prescription for characterizing the controlled evolutions in terms of time-independent 
effective Hamiltonians that would result in the same unitary propagator if applied over the same evolution interval. 
Imagine that a single control period T c consists of a cyclic sequence of \JC\ b.b. pulses, V — {Pk, Tfe}ji.J 1 , with 
n!S Pk = I- Then 

M M 

U{T C ) = exp(-iff e// T c ) = J] UlU {At m )U m = J[ exp(-iH m T m T c ) , (1) 

m— m— 

where the first equality defines the effective average Hamiltonian, and the "toggling-frame" transformed Hamiltonians 
H m are determined by the composite rotations U m = Iil'li-Pfcj ^ = — 1, Uq = 1. Here, M — \JC\ or 

M = \K\ + 1 depending on whether the sequence is arranged so as to allow evolution in the 1-frame in a single or 
a pair of control sub-intervals - being, in any case, U\jc\ = 1 by cyclicity. Because Pk G Gb.b., then G C Gb.b. as 
anticipated. However, G need not itself be a subgroup, neither does G = Gb.b. in general. For instance, allowing 
for G 7^ Gb.b. may be crucial for retaining universal control over decoupled dynamics Although the effective 
Hamiltonian ([!]) can be systematically calculated as a power series in the controllable parameter T c (Magnus series 
H ) , AHT is practically useful in the limit of fast control T c —> |^^,|| , where lowest-order contributions in T c suffice 
for an accurate description. In this limit, H e ff approaches 



H eff ^H=Y, nUlHUk = A g (H) , (2) 
fce/c 

with leading corrections accounted for by 

H 1 = ~ Wh ^ Hn ^ TmTn ' ^ 

c m>n 

Under the conditions ensuring convergence of the series defining H e ff, this correction is at least 0(T C ). The zero-th 
order approximation (^|) is, of course, exact if the various transformed Hamiltonians commute - in which case any 
time scale constraint disappears. 

In passing, it is worth noting that the most general transformation that AHT allows involves, as given in (j^), 
a weighted (convex) mixture of unitary operators. Extended to the space of linear operators End(Hs) over Hs, 
the action (|^) defines a trace-preserving, unital, completely positive map Ag. Examples abound in NMR where 
decoupling is achieved by control actions of this form. A representative case is the so-called whh-4 sequence used 
for homo-nuclear dipolar decoupling |Q. This corresponds to a sequence V = {P x , P- y , P y , P- x } of b.b. 7r/2-pulses, 
P a = cxp(— i7T(T a /4), and to = n = T3 = T4 = 1/6, T2 = 1/3 - which effectively averages out two-spin interactions 
proportional to 3 a\ <j{ — a 1 ■ over a cycle. 

The simplest realization of (||) occurs under the additional assumptions that control operations are equally separated 
in time and the composite rotations in G close a group (decoupling group ||). By letting Tk = 1/\G\, the quantum 
operation Ag of (||) reduces in this case to the projector on the centralizer Z(G) of G in End(7is), %{G) = {X £ 
End(H s )\[X,U k }=QVk}§§; 

11 fce/c 
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The resulting effective Hamiltonian acquires now a direct symmetry characterization as \H, Q] = - meaning that the 
controlled dynamics is symmetrized according to Q 0) , and the desired averaging action can be thought of as filtering 
out the dynamics that is not invariant under Q. 

One of the most famous (and practically useful) examples of a decoupling sequence with the above form is the 
Carr-Purcell (cp) sequence |l7]J|||, which in its basic variant is used for suppressing undesired phase evolution due to 
<j z terms - representing, for instance, applied field inhomogeneity in NMR. The sequence consists of repeated b.b. 7r a - 
pulses, a = x or y, separated by At, corresponding to a decoupling group Gcp a = {^,^a}, with 7r a = exp(— i7rcr a /2). 
Within the cycle time T c = 2At, the pulses can be arranged so as either t% = T2 = 1/2 (with M = \JC\ = 2) or 
T\ = T3 = 1/4, T2 = 1/2 (with M = |/C| + 1 = 3). Although the two sequences clearly lead to the same H as resulting 
from (||), the second option is actually superior in terms of the overall averaging accuracy - as it corresponds to a 

time-symmetric cycle for which all odd-order corrections H 2 +1 , £ = 0, 1, . . ., vanish 

In a closed-system setting, H$e — and the control problem consists in turning off (selectively or not) undesired 
contributions to Hs- Note that this amounts to simulating the reachable AHT Hamiltonians by H$. In particular, 
b.b. control may enable universal Hamiltonian simulation if an arbitrary effective evolution can be simulated by Hs- 
In an open-system setting, the focus is on averaging out the error generators {S a }, i.e. ensuring that Ag(J\f) = 0, in 
order to achieve decoherence suppression. The limit of fast control may be especially stringent to met in this case, as 
it requires that T c < r c for the shortest correlation time associated with the environmental noise |5|||[l0| - which can 
be prohibitively small in Markovian or quasi-Markovian dynamics. In both situations, the attainable control goals 
are influenced by two main factors: the available knowledge of the interactions to be manipulated and/or averaged 
out, and the overall available control resources. 

Some notable results about decoupling are worth mentioning. 

• Selective averaging (or refocusing) requires knowledge of the transformation properties of the interactions to be 
turned off with respect to the achievable decoupling sets Q . For a decoupling group Q , a necessary condition for 
selectivity is that the centralizer is non-trivial, Z(Q) ^ {Al} (A complex) [|]]|]. A variety of efficient schemes (i.e., 
polynomial in n) exist for selective decoupling in systems with at most bilinear interactions jl8| . These are applicable 
to NMR and NMR- like Hamiltonians in both the weak and strong coupling limit. 

• For arbitrary interactions (Hs or A/"), maximal averaging (or complete decoupling, or annihilation) is possible in 
principle by letting Q to be a unitary error basis on Hs ]q,pT|. Although this choice can be shown to be optimal, 
the resulting scheme is inefficient as the number of required b.b. control operations is \Q\ = N 2 = 4™ which grows 
exponentially with n. Again, efficient schemes can be designed if no terms higher than bilinear ones are known to be 
relevant |l8) . 

• Universal simulation of an arbitrary effective Hamiltonian on Hs can be implemented in various ways, depending 
on the set of actual simulation requirements and control resources. Suppose that, for a given Hs, a decoupling set 
Q exists, such that [Hs, Ag(Hs)] ^ 0, and that periods of free evolution under Hs can be alternated with periods of 
controlled evolution under Ag(Hs). Then, similarly to the twisted decoupler schemes discussed in any Hamiltonian 
L belonging to the Lie algebra generated by iHs, iAg(Hs) can be reached in principle - implying universal control in 
the generic case JTs| ] . Even in the unfavorable situation where Hs may consist of a single term (say a z for a qubit), 
arbitrary Hamiltonians can still be engineered if, for instance, both a set of b.b. operations averaging Hs, and a 
slow application of a Hamiltonian Z G Z(Q) can be effected. Then one can again alternate evolutions under Hs with 
controlled evolutions where Z is applied in parallel with the decoupler Q |8), obtaining as above universality in the 
generic case. (In the qubit example, with Hs = o z , one can choose Q = Gcp w and apply Z = a x via weak/slow control 
[|S|). Of course, such programming procedures may require additional external capabilities beyond b.b. control 
and, in general, they will not ensure universality starting from an arbitrary (possibly unknown) Hs- An elegant 
approach applicable to this general situation has been recently elaborated in pi| , where the possibility of arbitrary 
universal simulation is related to the identification of special decoupling groups called transformer groups. Starting 
from any Hamiltonian A, decoupling according to a transformer group is able to map A into any desired effective 
Hamiltonian. For instance, a transformer group for a two-dimensional system (a single qubit) is generated under the 
natural representation by the four b.b. operations {icr x , ia y , ia z , R}, R being the rotation by 27r/3 about the axis 
fi = (1, 1, l)/\/3 that cyclically permutes the Pauli matrices. 

III. ENCODED DYNAMICAL DECOUPLING 

So far, the full N = 2"-dimensional state space of n physical qubits has been exploited. However, restricting to 
a A^-dimensional quantum code, Nl < N, may prove extremely useful in QIP - the resulting benefits sometimes 
largely compensating the overheads and complications arising from dealing with a smaller number nj, < n of encoded 
qubits. In particular, the two primary motivations for seeking appropriate encodings are either to ensure protection 
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against noise - via active error-correcting codes |2Cj or passive noiseless codes on decoherence-free subspaces/noiseless 
subsystems (DFSs/NSs) |2l|- ^ , |l5| - or to allow for alternative routes to universality based on the physically available 
interactions on appropriately defined subsystem qubits - in the so-called approach of encoded universality |23|,|4| . 

A code Hl can be generally thought of as a distinguished subsystem of the physical state space of 5, determined 
by a correspondence of the form 

H s ^H L ®Hz®K, (5) 

for some Hz,TZ. Hl is the logical (or computational) factor, which reduces to a proper subspace Hl C Hs when the 
"syndrome" co-factor is one-dimensional, Hz — C. The summand 1Z collects the non-computational states in Hs- 
A code can be algebraically characterized with respect to a suitable algebra A of operators on Hs- The prototype 
example is the a noiseless code, whereby the appropriate algebra A is the associative interaction algebra |^2j containing 
the complex linear combinations of arbitrary products of Hs, all (or a subset of, see Sect. V) the error generators 
S a s, and the 1. Then in general A can be expressed, with respect to an appropriate basis in Hs, as a direct sum of 
dj-dimensional complex matrices, each appearing with a multiplicity nj, 

UAU^ = 01 nj ® Mat dj (C) , (6) 

JEJ 

where the change of basis U in Hs is made explicit, and ^2 Je j njdj = N. With respect to the same basis, the 
algebra A' = {X <G End(Hs) | [X, A] = 0} (commutant of A in End(Hs) H) represents as 

UA'U^ = 0Mat nj (C)(g)l dj . (7) 
JeJ 

Thus, under the unitary map U, the state space Hs becomes isomorphic to 

H s * Cj ® Vj ~ C nj ® C dj , (8) 

JeJ JeJ 

i.e. a direct sum of effectively bi-partite subspaces. Each of the left (or right) factors in this decomposition can 
be associated with a „4-code (or „4'-code), respectively. Typically, to obtain a noiseless code (a DFS or a NS), one 
selects a fixed factor Cj t = Hl, with dimension Nl = nj t - in which case, with respect to (||), Hz = T^M, and 
1Z = @j-Lj t Cj <g> T>j. Note that, because of ((?]), the code is an irreducible subspace of the commutant A' J^|J. 

If Nl > 2" L , then Hl can protect the state space of ul logical qubits against noise in A. While the global 
structure of Hl is sufficient for establishing storage or even existential universality results over Hl, an additional, 
crucial requirement on the local structure of Hl stems from the tensor product nature of QIP. In other words, an 
encoded tensor product structure on Hl is necessary for addressing notions of efficient simulation or universality 
over Hl- This issue is especially evident in the encoded universality approach, where the primary algebraic structure 
one considers is the Lie algebra £ generated under commutation by the set of easily implementable Hamiltonians 
p4j . (Formally, C plays the same role as A' in the noiseless approach.) Even if the code size Nl is large enough to 
accommodate many qubits, without a precise mapping that define encoded qubits there is no way for assessing the 
potential of this set of interactions in terms of one- and two-qubit encoded gates useful for implementing a quantum 
circuit. In its essence, properly defining this encoded tensor product structure is equivalent to properly constructing 
qubits in a given physical system |2j|. While no conclusive solution seems available to date, a practically motivated 
approach consists in identifying single-qubits encodings into small blocks of (2 to 4) physical qubits, and then inducing 
a tensor product structure by adjoining (or "conjoining") blocks |23|-|25f| . To do so, the system is partitioned into 
clusters {a} of physical qubits, i.e. Hs = JJ e H^ Ce \ and a mapping of the form (||) defines the state space Hl { of the 
£th encoded qubit starting from H^ CI \ An overall structure as in (||) still emerges with 

H L =H Ll ® ...» H LnL , (9) 

and TZ grouping all contributions involving TZi for at least one cluster. The code Hl can still be algebraically 
characterized as being, in general, embedded into one (or more) invariant subspaces of A 1 - or C, as appropriate 
pl] . p5|] . For definiteness, I focus on the case where encoding is motivated by noise protection against A henceforth. 
Let the set of encoded qubits be specified by encoded qubit observables, 

Wa e } , a = x,y,z; £ = l,...,n L , (10) 

satisfying Pauli-matrices commutation and anti-commutation rules p3|, and belonging to A'. The relevant situation 
for introducing encoded dynamical control assumes that the natural system dynamics - specified by Hs and possibly 
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by some error generators in AT whose effect is not eliminated by the encoding - are expressible in terms of encoded 
qubit observables. The goal is then to actively turn on/off selected encoded interactions without spoiling the benefits 
associated with the underlying encoding. Let U(Hl) — U(2 nL ) denote the group of unitary operators over Hl- As in 
the un-encoded case, the decoupling problem can be defined by a discrete set Q L — {U^} of b.b. control operations, 
which represent encoded rotations over Hl, Q l C U(Hl)- In the simplest setting, Q L will itself form a group of 
encoded rotations {encoded decoupling group). 

There are two minimal requirement for an operator Ufc to provide a legitimate unitary transformation over Hl'- 

• the gate should never draw the system outside the protected region; 

• the qubit mapping should be preserved at the end of the gate (2(| . 

Both conditions are satisfied if U% is generated by a Hamiltonian that is expressible in terms of encoded qubit 
observables - as it suffices to the present purpose. Accordingly, = exp(ii]SA^), for effective strength and time 
parameters 77, 8, respectively - in such a way that the limit 77 — > 00, 8 — » with a finite b.b. action r/S can be 
achieved. While specifying how rotations are effected is irrelevant in the idealized b.b. limit of instantaneous control 
actions considered so far, it clearly becomes important in a realistic scenario where pulses have a finite duration and 
the evolution during the pulses should therefore be taken into account. Within AHT, compensation schemes have 
been developed for dealing with pulse-length corrections [Q. For application with encoded pulses, it is necessary 
to preliminarily make sure that the benefits of the encoding are not lost during the pulses. Suppose that unitary 
operations in U(Hs) exist, whose action is not generated by Hamiltonians in A', but whose net effect matches, 
upon restriction to Hl, the one associated with some Ufc. Then besides the correction effects that also appear for 
exp(ir]8A^) for finite 8 (and A\ 6 A'), additional errors are associated with the departure from Hl during 8. While 
elimination of these effects motivates in principle the necessity of using encoded Hamiltonians, in practice different 
compensation techniques may be attempted, for instance by resorting to robust control design E ] . 

With these definitions and caveats in mind, the applicability of decoupling methods directly carries over to the 
encoded case, once qubits and qubit operators are formally replaced with their encoded counterparts. Thus, if an 
operator X = -F[{o^}] has a given structure in terms of physical qubit observables, and a decoupling scheme according 
to Q — {Uk} accomplishes a desired averaging effect via 

A s (X) = J2rkUlxU k , (11) 
keic 

then an equivalent averaging effect is obtained on an encoded operator X L with the same functional dependence 
X L = F[{<Ja e }] on encoded qubit observables, via the encoded quantum operation 

A g 4X L ) = Y,TkU^X L Ut. (12) 

keic 



IV. ENGINEERING OF ENCODED CLOSED-SYSTEM DYNAMICS 



Suppose that the noise generated by operators in Af is fully taken care of by the chosen encoding, or that noise 
is not a concern to begin with, as in the encoded universality approach. Then encoded dynamical decoupling may 
provide a tool for encoded universal Hamiltonian simulation. 

If, as assumed above, H$ expresses in terms of encoded observables, then the natural dynamics already implements 
a non-trivial logical transformation over the code. This provides the primary input to be exploited for engineering 
desired effective evolutions by encoded manipulations. As in the un-encoded case, detailed knowledge on the structure 
of H$ may or may not be a data of the problem. For arbitrary Hs, the results established in JlTJ ] imply that universal 
encoded simulation is achievable, in principle, if a finite encoded transformer can be constructed. The transformer 
groups so far identified |ll|l could be useful in principle for code size Nl = 2, 3 - although, unfortunately, practical 
impact is limited by the large number of encoded rotations involved (\G L \ — 24 for the above-mentioned single-qubit 
transformer, and \Q L \ > 168 for dimension 3). 

For the less ambitious target of generating a universal set of encoded evolutions starting from a given (known) 
Hs, simpler encoded programming strategies along the lines sketched in Sect. II may suffice. For instance, in the 
generic case, the Hamiltonian Hs, together with a non-commuting Hamiltonian Aql (Hs) obtained from Hs via some 
encoded decoupling procedure, will fulfill the conditions for generating the whole Lie algebra u(Hl) of anti-Hermitian 
encoded Hamiltonians - thereby implying universality over Hl, at least at the existential level. The required group 
Q L of encoded b.b. rotations may be as simple as an encoded Carr-Purcell-type group. Let, for instance, ir L denote an 
encoded it rotation (acting on one or more qubits), and Q^ P the associated encoded group. Then it is always possible 
to separate terms in Hs which are symmetric (s) and anti-symmetric (a) under Q^ p , 
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(13) 



with 

n L H i s s) 7r L = H s s = U gL (Hs), n L tt l = -H% . (14) 

Thus, the above argument generally applies provided [Hs, Ilgt (Hs)} = [H s ,Hg] ^ 0, and the two Hamiltonians 
Hs, Hg CI ,(Hs) accordingly suffice for universal encoded control. An experimental demonstration using a DFS-qubit 
encoded into the zero-quantum subspace of two nuclear spins is reported in E3 . 



A. Example: A single NS-encoded qubit 

A similar procedure could be relevant for obtaining universal encoded control over a NS-encoded qubit of three spin 
1/2 under general collective noise p^ |. Suppose that the noise generators S a are the global, permutation-invariant 
Pauli operators S a — X2j=i a a> a = x jVi z i an d that the natural system Hamiltonian has the isotropically-coupled 
form 

H s = ClS z + J12S12 + J23S23 + J31S31 , (15) 

where the Heisenberg exchange coupling Sjk — $ • <J k , and the parameters 0, Jjk are real. Then A is the algebra of 
completely symmetric operators over Hs — (C 2 )® 3 , and A' can be identified with the group algebra of the permutation 
group 53 (under the natural representation in Hs)- A NS code under A is identified by a correspondence of the form (H) 
- with Hz — C 2 carrying the irreducible representation of su(2) corresponding to total angular momentum J = 1/2, 
TLl — C 2 carrying the two-dimensional irreducible representation of £3, and 1Z ~ C 4 being the invariant subspace of 
states with total angular momentum J = 3/2 |^2|,^3|,|8) . Explicit expressions for encoded qubit observables (and the 
associated logical states) are given in [||). In particular, 

°*=l\ (l + aia), °y =l -^(s 23 -S3i) ■ (16) 

This allows to rewrite the Hamiltonian ( [To] ) , up to irrelevant contributions which are constant over Hl , as 

H s =l (2 J12 - J23 - J31) ° L X + \/3 ( J31 - J23) Py ■ (17) 

Note that the vanishing of Hs for a fully symmetric coupling network, J.^ = J Vj, k, correctly verifies the identity 
action of permutation- invariant operators over TCl (Hs G A in this case). Given Hs in the form (|l7|), the above- 
mentioned universality scheme based on encoded Carr-Purcell sequences becomes applicable in principle provided 
one has to ability to enact rapid, encoded n L pulses. Because, with respect to the chosen code, the action of 
the operator is effectively identical to swapping the physical qubits 1 and 2, this is for instance achievable if 
the exchange Hamiltonian S12 can be switched on for the appropriate time. As a result, one generates an effective 
encoded Hamiltonian n g i (Hs) = (2J12 — J23 — Jsi)^ which, together with Hs, allows for universality. Although not 
directly applicable to the weakly-coupled molecule used in |27j to realize the above NS, and generally demanding for 
NMR implementations as the J-coupling parameters are not controllable, these ideas could prove viable for solid-state 
implementations where Heisenberg interactions are in principle fully tunable p4}| . Explicit universality constructions 
relevant to quantum computing architectures based on exchange interactions are provided in Jl^| . 



B. Example: A pair of DFS-encoded qubits 

As a further example, motivated from NMR, imagine two pairs of spin 1/2 nuclei, corresponding to different species 
(Hydrogen and Carbon, for instance), subjected to block-collective dephasing noise. Let the two clusters be associated 

with the pairs (1, 2) and (3, 4), respectively. Then the two error generators are = a\ + a^, — c 3 + er 4 , and 
A = A^z^ ®A^z 2 ^ in terms of the interaction algebras A z c ^ generated by S z Ce "* for dephasing on 2 qubits H] . Protection 
against errors in A can be accomplished by replacing each physical pair with a DFS-encoded qubit supported by the 
states having — and S Z C2 ^ — i.e., 

|»L I J £a > = |i£ 1 >®|iO, IOl,) = |01) (C{) , |l ii ) = |10)^), £=1,2. (18) 
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In this case, with reference to the general structure (||), each DFS encoding has the form ~ H.L e © Hi, with 

7i.Lt = s P an {|0-Lf)i — C 2 , TZt — span{|00)£, |11)^} — C 2 (Hz e — C is irrelevant for subspace-encodings) . 

Correspondingly, for 4 spins Hl — C 2 ® C 2 is a 4-dimensional subspace of the 6-dimensional zero-quantum subspace 
corresponding to total z angular momentum S z — Y) j a J z = 0, j = 1, . . . , 4, and 1Z collects all contributions where 

S { z Ci) ^ for at least one pair. Encoded observables for the above qubits are provided, for example, by the choice llj] 



^ =L 2 



\ ft - ^) , ^ =L \ (o\al + alafj , (19) 



and similarly for qubit L2. 

In general, the physical spin system will exhibit a strongly-coupled spectrum described by an internal Hamiltonian 
composed of both Zeeman spin-field and indirect spin-spin interactions, 

H S = + E \^'^ • Y.1'1 ■ ( 20 ) 

j'=l,— A j<j'=i,...,i j j<j' 

where the chemical shifts and J-coupling parameters are understood in frequency units. For usual values of the 
static Zeeman field, the contribution to the total energy of a given pair of spins due to the coupling Hj., 

can be treated as a perturbation with respect to Tfz ., - the diagonal o\o\ and off-diagonal + cr^ y a 3 y terms 

leading to first- and second-order correction effects in Jjj> , respectively. Because the differences in the chemical shifts 
\vj — Vji I are larger when different nuclear species are involved, the approximation of neglecting off-diagonal couplings 
(weak-coupling limit) is well justified for hetero-nuclear interactions. Thus, (^0|) can be effectively replaced by 

H s = ^ j ai + ^(j 12 a 1 -a 2 + J 3i a 3 -a i ) + Y \ J ji<°i- (21) 

. I I j=l,2;j'=3,4 



In terms of the encoded qubit observables given in (19), the chemical shift terms immediately rewrite as combinations 
of S z c ^ operators (constant over Hl) and logical a z e operators, whereas the homo-nuclear spin-spin interactions 
contribute with logical of f o per ators (and additional constant terms proportional to o^o^). The remaining hetero- 
nuclear bilinear couplings in ( |2l| ) give 



with coefficients 

A 
C 



ASi Cl) Si C2) + BS^a% 2 +CSi C2) a^ + Da^a^ 2 , (22) 

g (j%3 + Jl4 + -hi + ^24^ , B = — ^ J13 — J14 + J23 — ^24^ , 

J (jl3 + Jl4 — J23 ~ ^24^ , D = — ^Jl3 — J14 — J23 + -^24^ • (23) 



Again, an overall identity action over TLl is obtained for a symmetric coupling network - as well as for other special 
coupling patterns in the system. Assuming that none of these non-generic circumstances is met, the action of the 
overall Hamiltonian Hs over the code is finally: 



Ha 



it (Ai/12 of 1 + Ai/34 of 2 + J 12 of 1 + J 34 of 2 + D of 'of 2 ) , (24) 



with Az/12 = v\ — v%, A^34 = 1/3 — 1/4, respectively. 

Remarkably, the structure of the Hamiltonian (|24|) for the two encoded qubits is very similar to the one of the 
Hamiltonian describing two weakly-coupled physical spins - except, in the ordinary NMR setting, control along 
the transverse a x (or a y ) directions is directly supplied by external radio- frequency fields. While the potential of 
( p4| ) for universal encoded control should be expected by analogy with the un-encoded case, various schemes may 
be specified depending on the actual control capabilities. At the existential level, universality follows from the 
ability of rapidly effecting a single encoded rotation, say a "hard" (non-selective) encoded ir L pulse, 7if = Tif x 7!f 2 
about the encoded x axis. Because the associated averaging produces the encoded effective Hamiltonian U s l (Hs) = 
J12&X 1 + J34&X 2 + Da z 1 a z 2 , by appropriately alternating evolution periods under Hs and Hql (Hs), arbitrary encoded 
Hamiltonians can be enacted in principle through repeated commutation p9[ . 

Constructive prescriptions become possible as soon as a wider range of control options is accessible. In particular, 
schemes for effecting a universal set of encoded gates - including (i) all single qubit rotations and (ii) a two-qubit 
phase coupling - can be provided. Assume that, in formal analogy with the standard NMR setting, the parameters 
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J12 and J34 are independently tunable at will. To selectively rotate each encoded qubit about the x axis, turn on 
the appropriate term in and simultaneously refocus a^ e evolutions by applying encoded n^ir^ 2 pulses as 
above. Note that although this action preserves the erf 1 erf 2 coupling, the corresponding evolution can be neglected 
under the usual assumption that single-qubit rotations can be effected rapidly enough. In fact, the whole evolution 
induced by ( |24| ) can be neglected during the x rotation if the strength of the appropriate J12 or J34 parameter can be 
arbitrarily controlled as assumed. Also, note that the same sequence of hard ir^ 1 ^ 2 pulses allows to selectively leave 
this phase coupling alone if both J12 and J34 are kept to zero - except during the (infinitesimal) durations of the 
refocusing pulses. Selective rotations about the appropriate encoded z axis are only slightly more demanding, as they 
require the hard logical ir x pulses to be replaced by "soft" encoded i r^ e on the intended spin alone. The relevance of 
these and related universality constructions has been highlighted in . 

However, the above tunability requirements are too stringent, for instance, for NMR implementations - where 
the natural Hamiltonian ( pi] ) is always active, and interactions can be only effectively set to zero via appropriate 
refocusing. Additional constraints may arise for typical implementation parameters. In particular, one may expect 
that the encoded chemical shift evolutions will dominate in (p4f), i.e. A^i2,A^34 ^> J12, J34 3> D. Nevertheless, 
universal control can still be gained by having access to a small set of encoded operations. No ability of arbitrarily 
tuning encoded Hamiltonians is assumed, apart from the fixed combinations of strength and times corresponding to 
logical ir L pulses. However, unlike the previous case, encoded n L rotations about two axes are generally needed. The 
encoded decoupling schemes for implementing the required universal set of coupling are as follows. Selective rotations 
about the appropriate z axis are easiest, as one can take advantage of the natural averaging of the <r x e and a^a x 2 
that is enforced by the above hierarchy. One need only refocus the undesired encoded phase evolution (say, cf 2 ) 
by applying sequences of encoded selective n x e pulses (say, n x 2 ). It is easily seen that decoupling purely based on 
simultaneous tt^tt^ 2 pulses as above is not useful, without the ability to separately controlling the coefficients in the 
resulting effective Hamiltonian, to implement selective x rotations or the two-qubit coupling. However, a selective 
implementation of, say, erf 1 can be engineered, for instance, by the encoded pulse sequence 

tt^tt^ - At - - At- nt^t - At - tt^tt^ , At = T c /4 , (25) 

which corresponds to subjecting qubit 2 to maximal encoded averaging according to Q^ nax = erf, er^ , erf}, and 
cycling twice qubit 1 through the encoded Carr-Purcell sequence refocusing erf 1 . A similar sequence works for rotating 
the second encoded qubit, by interchanging 1 and 2 in (H|). If all single-qubit interactions are refocused instead, by 
applying the above Gm ax to both qubits, 

TT^TT^ - At - TT^Trf 2 - At- 7#7# - At - Tr^Tr^ , At = T c /4 , (26) 

the erf 1 erf 2 evolution is selectively extracted from (|24|). 

The problem of how to actually effect the required logical ir L rotations remains challenging in reality, as the 
evolutions induced by the available radio- frequency fields do not in general correspond to encoded Hamiltonians. 
Thus, even though it is always possible to enforce a unitary propagator whose net action is the same as the one of 
the required encoded rotation, the corresponding sequence of physical gates may be complicated and, as noted above, 
will cause departure from the code during the finite pulsing time. A possible way out is explored in JL4| , based on the 
idea of compensating for the resulting exposure to noise by optimizing the length of the relevant control sequences 
and by incorporating intrinsic robustness features using composite pulse techniques p9[ . The situation is relatively 
straightforward for the special case of ir L pulses, as the encoded rotations they effect have a simple translation in 
terms of realizable physical spin propagators. The simplest instance is the action of a Trf^f 2 , which is identical to 
the action induced over Hl by a hard n x pulse on all four spins. In fact, the following correspondence for the logical 
w L pulses involved in the above encoded sequences holds (up to irrelevant phase factors): 
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Accordingly, universal manipulation of encoded evolutions with reduced error rate is still achievable if robust ways 
for effecting the above physical operations can be devised |Q|o) ■ 
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V. ENGINEERING OF ENCODED OPEN-SYSTEM DYNAMICS 



A further direction for application is suggested by looking at encoded dynamical decoupling as a tool for encoded 
quantum error suppression. Similar to the idea of concatenating passive noise protection via DFSs/NSs with quantum 
error correction - by operating error-correcting codes directly on DFS/NS-encoded qubits J3l]], the idea is now to 
concatenate passive noise protection with active dynamical control - by effecting error suppression directly on encoded 
qubits. 



A. Example: A decohering encoded qubit 

By analogy to the prototype example of a single decohering qubit analyzed in p] , the prototype situation is provided 
by a single decohering encoded qubit. Imagine a system S composed of two physical qubits, which are subjected to a 
purely dephasing interaction as a result of the coupling to two bosonic environments B\, B2 i.e., the overall system 
is described by a Hamiltonian of the form 

H = H s + H Bl + H B2 + H SBl + H S b 2 , (27) 

for uncoupled subsystem's Hamiltonians 

Hs=J2 Y ai ' Hb > + ^ = £ + E -4 2) ^ 2)t , (28) 

3=1,2 k i 

and a total interaction Hamiltonian 

h SBi + H S B 2 = ^{EW+fi+EW+f))- ( 2g ) 
3=1,2 k i 

In the above equations, u>j, j = 1,2, are the natural frequencies of single-qubit evolutions, ui >&jt are bosonic 
operators for the environment B\ (similarly for B2), and the parameters g^ (de^) determine the coupling strength 
of qubit j to mode k of bath 1 (or mode £ of bath 2, respectively). Identity operators on the appropriate subsystems 
are also understood as needed. Then implementing decoupling according to the group Gcp x — {^j&x&x} would 
dynamically suppress both error generators - thereby allowing to reduce the dephasing noise on both qubits to a 
level, in principle, as low as desired. However, as recalled earlier, this is only effective in practice in the limit of fast 
control where pHlPf 

r^min^rW}, (30) 

Tc denoting the (shortest) correlation time of the ith environment. Suppose now that B2 is sufficiently "slow" , but B\ 

is "fast" , so that decoupling at the rate determined by becomes unfeasible. Although implementing decoupling at 

the slower rate set by B2 reduces the error rate originating from B2, there is no actual guarantee that the overall error 

(2) 

rate is suppressed due to the possibility of decoherence acceleration from modes at frequency higher than ~ 1/tc 
in B\ [5j. In any event, both qubits would remain effectively exposed to noise. Are there other noise control options 
worth being considered? 

For arbitrary dephasing and just two qubits, quantum error correction does not help, neither does passive noise 
protection - unless some symmetries can be identified in the noise. Suppose that, in the above interaction Hamiltonian 
(p9|), gp = g^ — to a good accuracy. Then the overall coupling rewrites as 

H SBl + H SB2 = S z <g> + al ® B[ 2) + a 2 z ® B { 2 2) , (31) 

for S z = al + a 2 z as earlier defined, and appropriate environment operators = Y^kSkfik ^fc^)' ~ 

Y^tQij {"e + "i )i 3 — 1, 2. Physically, this situation corresponds to a noise process consisting of both "fast" 
collective dephasing, due to B\, and "slow" independent dephasing, due to £?2- Whenever a symmetry in the open- 
system dynamics is present, major gains should be expected by seeking for an appropriate encoding into a DFS or 
a NS. In our case, the relevant DFS is the one already introduced in Sect. IVB i.e., the one spanned by states 
corresponding to total zero angular momentum along z, S z — 0. Thus, a logical DFS qubit is defined by the encoding 
given in (f||) (for just one qubit), and 







Hl = span{|0 L ), |U)} = span{|01), |10>> 



(32) 



By construction, this encoded qubit is perfectly protected (with infinite distance) against the noise due to B\. However, 
mixing with the degrees of freedom of J3 2 is still induced through the error generators , j = 1 , 2 in ( |3l|) . The key 
observation is that this residual noise preserves the coding space, corresponding to a purely decohering coupling 
between the encoded qubit and B2. Because, by using (|l9|), the error generators simply express in terms of encoded 
observables 

a z ~L Y +a * ' 2 =L ~~2 z ' ' > 

the action of the total Hamiltonian H ( p7j ) on Hl finally rewrites as 

H= L Aua^ + H B +(J L z ®B Zl (34) 

f2) (2) 

where Aw = [oji — W2)/2, Hb = Hb x + Hb 2 , and B z = [B\ — B 2 )/2- This form makes it clear that the action of H 
on the encoded qubit is formally identical to the action of the diagonal spin-boson Hamiltonian on the physical qubit 
considered in flp^|. 

By the same argument valid for suppressing decohcrcncc on a single physical qubit, implementing a sequence 
of equally spaced, encoded pulses sufficiently fast with respect to the slower rate determined by is now 

guaranteed to suppress the encoded error rate by a factor of (at least) 0[{T c /tc 2 ^) 2 ] ||. This corresponds to encoded 
dynamical decoupling according to Gcp x = 17 x }j where a x is again given in (p9|). Thus, if the required encoded 
71"^ rotation can be effected, suppression of the "slow" dephasing from B2 can be accomplished without re-introducing 
exposure of the encoded qubit to the "fast" dephasing from B\. 



B. Generalizations 



Some generalizations of the above example are worth pointing out. First, the same encoded decoupling scheme is 
effective at suppressing encoded phase errors in a situation where the internal two-qubit Hamiltonian is governed by 
a generally anisotropic exchange Hamiltonian 

Hs = + Y a " + J &* a * + J »y + J ™ a * a * ' (35) 

with Jf 2 = Jyj, = J anc ^ either Jf 2 arbitrary, Jf 2 = 0, or Jf 2 = J (the so-called XXZ, XY, and isotropic models, 
respectively ]13|] ). In the isotropic case in particular, which was also examined in Sect. IV and is directly relevant to 
NMR QIP, the encoded open-system Hamiltonian ([54j) is modified to an encoded spin-boson Hamiltonian [^|| as 

H' = L Aua^ + Ja L x + H B + cr, ® B z . (36) 

Once encoded suppression of the a z error generator is accomplished as above, control according to the Hamiltonian 
can be re-introduced, if desired, by implementing one of the programming schemes described in Q - implying the 
possibility of retaining universal noise-suppressed encoded control. 

Situations involving hybrid dephasing processes with highly correlated components and slow residual noise in multi- 
qubit systems may be also relevant to NMR. If, for instance, in a four-qubit system, dephasing from B\ is pair-wise 
correlated on qubits (1,2) and (3,4) as considered in Sect. IVB, then encoding into the tensor product of the two 
DFSs described there ensures protection against B\. Once this is done, applying sequences of encoded tt^tt^ 2 pulses 
at the appropriate rate causes a suppression of any (independent or correlated) residual dephasing due to B^- 

Finally, similar ideas may be more generally applicable to situations where the encoded error generators for the 
residual noise are expressible through encoded observables. If, for instance, starting from the above two-qubit Hamil- 
tonian (^), the action over Hl may be cast into the form 

H = L Alu + H B + (T z ® B z + ® B x , (37) 

for appropriate operators B z _ x on environment B 2 , then the encoded qubit suffers from depolarizing noise. The result- 
ing error rate can be suppressed, in principle, by using the encoded annihilator Gmax = {l 1 ,^,^,^} introduced 
above. Note that, in terms of coupling to the physical degrees of freedom, the interaction ( ]37j ) involves multiple-qubit 
excitations as error generators. Interestingly, system-bath couplings allowing for similar multi-qubit processes to low- 
est order are the only class of interactions, in addition to the ones exhibiting spatial symmetry, for which DFSs are 
known to exist |Q. 
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To summarize, if noise has both slow and fast components, so that achieving full decoupling of the physical degrees 
of freedom becomes unfeasible, then encoded dynamical decoupling may be useful in situations where a dominant 
symmetry in the fast noise can be exploited to obtain encoded qubits, and the generators for the residual slow noise 
are expressible in terms of encoded observables. In this event, using encoded rather than physical degrees of freedom 
allows the system to benefit already from enacting decoupling operations at the slower rate. However, this looser 
constraint on control time scales competes, as already emphasized, with tighter symmetry constraints on the useful 
control Hamiltonians - as encoded Hamiltonians are demanded, and they are not always easily available in physical 
systems. 



VI. CONCLUSIONS 



I have analyzed the relevance and potential usefulness of active dynamical control, as inspired by NMR spectroscopy, 
in the light of the QlP-motivated notion of encoded degrees of freedom. Ultimately, the resulting approach of 
concatenating active control with encoded qubits naturally stems from the program of regarding the information- 
carrying subsystems as the primary degrees of freedom for realizing QIP. Applications of encoded dynamical decoupling 
to control both Hamiltonian and non-Hamiltonian encoded evolutions have been considered in principle. While the 
actual viability within both settings strongly depends on the details of the specific implementation and the available 
control resources, I hope that these ideas will motivate further investigation and serve as guiding principles to further 
expand our capabilities to manipulating quantum systems and quantum information. 
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